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Nanofluids  have  evoked  immense  interest  from  researchers  of  multi-disciplines  from  around  the  globe 
due  to  their  fascinating  thermophysical  properties  and  numerous  potential  benefits  and  applications  in 
important  fields  such  as  microelectronics,  microfluidics,  transportation,  and  biomedical.  However,  there 
are  many  controversies  and  inconsistencies  in  reported  arguments  and  experimental  results  on  various 
thermal  characteristics  such  as  effective  thermal  conductivity,  convective  heat  transfer  coefficient  and 
boiling  heat  transfer  rate  of  nanofluids.  As  of  today,  researchers  have  mostly  focused  on  anomalous  ther¬ 
mal  conductivity  of  nanofluids.  Although  investigations  on  boiling,  droplet  spreading,  and  convective 
heat  transfer  are  very  important  in  order  to  exploit  nanofluids  as  the  next  generation  coolants,  consid¬ 
erably  less  efforts  have  been  made  on  these  major  features  of  nanofluids.  In  this  paper,  these  important 
cooling  features— boiling,  spreading,  and  convective  heat  transfers  of  nanofluids  are  presented  together 
with  exhaustive  review  of  research  and  development  made  in  these  areas  of  nanofluids. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

With  an  ever-increasing  thermal  load  due  to  smaller  fea¬ 
tures  of  microelectronic  devices  and  more  power  output,  cooling 
for  maintaining  desirable  performance  and  durability  of  such 
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devices  is  one  of  the  most  important  technical  issues  in  many 
high-tech  industries.  The  conventional  method  to  increase  the 
cooling  rate  is  to  use  extended  heat  transfer  surfaces.  How¬ 
ever,  this  approach  requires  an  undesirable  increase  in  the  size 
of  the  thermal  management  system.  In  addition,  the  inherently 
poor  thermal  properties  of  traditional  heat  transfer  fluids  such 
as  water,  ethylene  glycol  (EG)  or  engine  oil  (EO)  greatly  limit 
the  cooling  performance.  Thus,  these  conventional  cooling  tech¬ 
niques  are  not  suitable  to  meet  the  cooling  demand  of  these 
high-tech  industries.  Although  thermal  conductivity  of  a  fluid  plays 


S.M.S.  Murshed  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  15  (201 1 )  2342-2354 


2343 


Nomenclature 

CSf  an  empirical  coefficient  in  Rohsenow  correlation 
cp  specific  heat  (J/kg  K) 

D  tube  diameter  (m) 

/  friction  factor 

g  acceleration  of  gravity  (m/s2 ) 

h  heat  transfer  coefficient  (W/m2  K) 

k  thermal  conductivity  (W/m  K) 

I<  a  constant  in  Zuber  correlation 

L  length  (m) 

m  mass  flow  rate  (kg/s) 

q  supplied  heat  (W) 

q "  heat  flux  (W/m2) 

T  temperature  (K  or  °C) 

Greek  symbols 

0  particle  volume  fraction 

p  density  (kg/m3) 

/> i  dynamic  viscosity  (kg/m  s) 

a  surface  tension  (N/m) 

Subscripts 
b  bulk 

/  fluid 

g  vapor  phase 

i  inner 

o  outer 

nf  nanofluid 

p  particle 

s  saturated  liquid 

v  laminar  sublayer 

w  wall  surface  (heater  or  tube) 

Nondimensional  numbers 
Nu  Nusselt  number 

Pe  Peclet  number 

Pr  Prandtl  number 

Re  Reynolds  number 

Abbreviations 

BHF  burnout  heat  flux 

BHTC  boiling  heat  transfer  coefficient 

CHF  critical  heat  flux 

CNT  carbon  nanotube 

DNB  departure  from  nucleate  boiling 

Ni-Cr  nickel-chromium  (nichrome) 

ONB  onset  of  nucleate  boiling 


a  vital  role  in  the  development  of  energy-efficient  heat  trans¬ 
fer  equipment  it  is  known  that  fluids  possess  order-of-magnitude 
smaller  thermal  conductivity  than  metallic  or  nonmetallic  particu¬ 
lates.  For  example,  thermal  conductivities  of  water  (0.607  W/m  I<) 
and  engine  oil  (0.145  W/m  K)  are  about  5000  times  and  21,000 
times,  respectively  smaller  than  that  of  carbon  nanotubes  (e.g. 
3000  W/m  K  for  multi-walled  carbon  nanotubes,  MWCNTs).  Fig.  1 
compares  thermal  conductivities  of  various  commonly  used  liq¬ 
uids  and  materials  at  room  temperature.  It  can  be  seen  from  Fig.  1 
that  most  of  the  commonly  used  metallic  and  nonmetallic  par¬ 
ticulates  posses  orders  magnitude  higher  thermal  conductivity 
than  those  of  liquids.  Fig.  1  also  shows  that  metals  have  higher 
thermal  conductivity  compared  to  ceramics  (or  metallic  oxides). 
Therefore,  the  thermal  conductivities  of  fluids  that  contain  sus- 
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Material 


Fig.  1.  Thermal  conductivities  of  commonly  used  liquids  and  materials  at  room 
temperature  (some  small  thermal  conductivity  values  are  shown  above  the  numbers 
of  the  liquid  or  material)  [1-6]. 


pended  metallic  or  nonmetallic  particles  or  tubes  are  expected 
to  be  significantly  higher  than  those  of  traditional  heat  transfer 
fluids. 

Although  nanoparticles-suspensions  were  used  in  heat  trans¬ 
fer  studies  as  early  as  1984  by  Yang  and  Maa  [7]  and  then 
in  1993  by  Masuda  et  al.  [3],  it  is  only  in  1995  that  Choi  [8] 
at  Argonne  National  Laboratory  of  USA  coined  the  concept  of 
“nanofluid”  to  meet  the  cooling  challenges  facing  many  high- 
tech  industries.  Apart  from  Masuda  et  al.  [3],  there  are  other 
researchers  who  also  used  the  term  “nanofluid”  or  were  involved 
with  research  on  thermal  conductivity  of  nanoparticle  suspen¬ 
sions  during  1993.  Gass  et  al.  [9]  from  Switzerland  used  the 
same  term  “nanofluid”  to  express  minute  volume  of  fluid  (nano¬ 
liter)  in  microfluidics  study  in  1993.  A  German  researcher,  Grimm 
[10]  also  won  a  German  patent  on  the  enhanced  thermal  con¬ 
ductivity  of  nano-  and  micro-sized  particles  suspensions  in  the 
same  year.  Aluminum  particles  of  80  nm  to  1  p>m  were  suspended 
into  a  fluid  and  about  100%  increase  in  the  thermal  conductiv¬ 
ity  of  the  fluid  for  loadings  of  0.5-1 0vol.%  was  reported  in  his 
patent. 

The  definition  of  this  new  class  of  heat  transfer  fluids,  i.e. 
nanofluids  is  the  suspensions  of  nanometer-sized  solid  parti¬ 
cles,  rods  or  tubes  in  traditional  heat  transfer  fluids.  Nanofluids 
were  found  to  exhibit  significantly  higher  thermophysical  prop¬ 
erties,  particularly  thermal  conductivity  and  thermal  diffusivity 
than  those  of  base  fluids  [11-16].  Nanofluids  have  attracted  great 
interest  from  the  research  community  due  to  their  enhanced  ther¬ 
mal  performance,  potential  benefits  and  applications  in  numerous 
important  fields  such  as  microelectronics,  microfluidics,  trans¬ 
portation,  manufacturing,  medical,  and  so  on.  Recent  record  shows 
(Fig.  2)  that  there  is  an  exponential  growth  of  annual  research 
publications  on  nanofluids.  According  to  ISI  Web  of  Knowledge 
searched  results  (Fig.  2),  total  153  nanofluids-related  publications 
have  already  appeared  in  the  first  half  of  year  2010  (till  July 
14).  It  is  also  believed  that  there  are  more  than  300  research 
groups  and  companies  worldwide  who  are  involved  with  nanoflu¬ 
ids  research  [17].  Compared  to  studies  on  thermal  conductivity, 
few  works  have  been  reported  on  boiling,  droplet  spreading,  and 
convective  heat  transfer  characteristics  of  nanofluids  although 
these  features  are  very  important  in  order  to  exploit  nanofluids 
as  the  next  generation  coolants.  It  is,  therefore  imperative  and 
timely  to  provide  a  state-of-the-art  review  on  the  advances  in  boil¬ 
ing,  droplet  spreading,  and  convective  heat  transfer  research  of 
nanofluids. 
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Publication  year 


Fig.  2.  Annual  publications  on  nanofluids  since  2000  (publications  including  all 
types  of  journal  and  conference  articles,  patent,  news,  letter  and  other). 
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2.  Boiling  and  droplet  spreading  of  nanofluids 

2.2.  Boiling  heat  transfer— basic 

Boiling  is  a  phase  change  (from  liquid  to  vapor)  process  in  which 
vapor  bubbles  are  formed  either  on  a  heated  surface  or  in  a  super¬ 
heated  liquid  layer  adjacent  to  the  heated  surface.  Although  boiling 
is  a  complex  and  elusive  process,  it  is  very  efficient  mode  of  heat 
transfer  in  various  energy  conversion  and  heat  exchange  systems  as 
well  as  cooling  of  high  energy  density  electronic  components.  There 
are  two  common  types  of  boiling:  pool  boiling  and  flow  or  forced 
convective  boiling.  Pool  boiling  refers  to  boiling  on  a  heated  surface 
(heater)  submerged  in  a  pool  of  initially  quiescent  liquid  while  flow 
boiling  is  boiling  in  a  flowing  stream  of  fluid,  where  the  heating  sur¬ 
face  may  be  the  channel  wall  confining  the  flow.  Heat  flux  in  any 
boiling  is  one  of  the  most  important  parameters  in  designing  and 
operating  the  heat  transfer  equipments  with  high  heat  flux  such  as 
boiler,  evaporator,  electronic  equipment  and  rocket  engines.  Since 
there  are  very  few  efforts  are  made  on  flow  boiling  of  nanofluids, 
this  mode  of  boiling  heat  transfer  will  not  be  elaborated  further 
instead  pool  boiling  heat  transfer  is  briefly  discussed  here. 

In  pool  boiling  problems,  the  fluid  is  initially  quiescent  near 
the  heating  surface,  and  subsequent  fluid  motion  arises  from  free 
convection  and  the  circulation  induced  by  bubble  growth  and 
detachment.  Two  main  parameters  in  boiling  problems  are  the 
degree  of  wall  superheating  or  excess  temperature  (the  difference 
between  the  wall  and  the  liquid  saturation  temperature  at  the  local 
pressure,  AT  =  TW  -  Ts)  and  heat  flux  ( q ").  The  classical  pool  boil¬ 
ing  curve,  which  is  as  a  plot  of  heat  flux  versus  excess  temperature 
(AT),  is  shown  in  Fig.  3.  As  the  value  of  AT  increases,  the  curve 


traverses  four  different  regimes  or  modes  of  pool  boiling  (Fig.  3): 
(I)  free  convection,  (II)  nucleate  boiling,  (III)  transition  boiling,  and 
(IV)  film  boiling. 

As  illustrated  in  Fig.  3,  up  to  point  A  is  the  free  convection  regime 
which  occurs  since  there  is  insufficient  vapor  to  cause  active  boil¬ 
ing.  Small  temperature  differences  exist  in  the  liquid,  and  heat  is 
removed  by  free  convection  to  the  free  surface.  At  point  A,  isolated 
bubbles  initially  appear  along  the  heating  surface  and  this  is  termed 
as  onset  of  nucleate  boiling  (ONB).  Nucleate  boiling  occurs  between 
points  A  and  C.  At  this  regime,  isolated  bubbles  appear,  and  heat  is 
transferred  mainly  from  the  surface  to  the  liquid.  As  AT  increases 
(B-C),  more  nucleation  sites  become  active  and  bubbles  coalesce, 
mix,  and  ascend  as  merged  jets  or  columns  of  vapor.  At  point  C, 
the  maximum  heat  flux  or  critical  heat  flux  (CHF)  occurs.  Between 
points  C  and  D,  transition  boiling  occurs.  An  unstable  (partial)  vapor 
film  forms  on  the  heating  surface,  and  conditions  oscillate  between 
nucleate  and  film  boiling.  Intermittent  vapor  formation  blocks  the 
liquid  (higher  conductivity)  from  contacting  the  surface  (lower¬ 
ing  the  surface  heat  flux).  Film  boiling  occurs  beyond  point  D.  In 
addition  to  conduction  and  convection,  heat  transfer  by  radiation 
is  important  at  these  high  wall  superheat  levels.  A  stable  vapor  film 
covers  the  surface  in  this  regime. 

There  are  numerous  factors  such  as  heater  or  channel  surface 
conditions  (smooth  or  rough),  heater  size,  shape,  material,  diameter 
and  orientation,  degree  of  surface  wetting,  subcooling,  inclusion  of 
surfactants,  and  properties  of  liquid  affecting  heat  transfer  perfor¬ 
mance  and  bubble  dynamics  of  pool  and  flow  boiling.  These  factors 
are  widely  studied  both  theoretically  and  experimentally  [18-23]. 
For  instance,  effects  of  several  parameters  such  as  tube  diame¬ 
ter,  length,  surface  roughness,  and  inclination  on  pool  boiling  heat 
transfer  have  been  investigated  systematically  by  Kang  [21  -23  ].  His 
results  showed  that  there  are  significant  individual  or  combined 
effects  of  these  factors  on  boiling  heat  transfer  characteristics. 

2.2.  Boiling  heat  transfer  of  nanofluids— literature  survey 

The  boiling  critical  or  burnout  heat  flux  enhancement  of 
particle-suspensions  depends  on  the  particle  concentration,  pH  of 
the  solution  as  well  as  on  the  deposition  of  the  particles  on  the  sus¬ 
pended  heater  surface.  It  is  long  back  proven  that  addition  of  solid 
particle  in  base  fluid  can  alter  its  boiling  heat  transfer  performance. 
For  example,  Yang  and  Maa  [7]  first  used  nano-sized  A1203  particles 
as  small  as  50  nm  in  water  to  study  the  pool  boiling  heat  transfer 
characteristics.  They  reported  significant  increase  in  pool  boiling 
performance  for  very  small  volumetric  concentrations  (0.1 -0.5%) 
of  nanoparticles.  After  nanofluids  emerged  a  growing  number  of 
research  groups  have  got  involved  with  boiling  heat  transfer  study 
of  nanofluids  and  it  is  timely  to  review  their  research  findings. 

For  Al203-water  nanofluid  with  a  flat  plate  heater,  You  et  al. 
[24]  observed  a  three-fold  increase  in  critical  heat  flux  (CHF).  For 
silica  nanofluids,  similar  three-fold  enhancement  in  CHF  was  later 
reported  by  Milanova  and  Kumar  [25].  Whereas,  Das  et  al.  [26] 
reported  deterioration  of  boiling  heat  transfer  of  water  in  the  pres¬ 
ence  of  A1203  nanoparticles  in  it.  Their  outcome  was  partially 
attributed  to  the  properties  of  the  nanofluid,  boiling  surface  and 
the  interaction  between  the  two.  In  contrast  to  Das  et  al.  [26],  Wen 
and  Ding  [27]  showed  that  the  enhancement  of  pool  boiling  heat 
transfer  of  the  same  7~A1203 /water-based  nanofluid  was  about  40% 
at  1.25  wt.%  of  particle  loading.  In  another  study,  Witharana  [28] 
investigated  the  boiling  heat  transfer  performance  of  two  types 
of  Au  and  Si02-laden  aqueous  nanofluids  in  a  cylindrical  vessel 
under  atmospheric  pressure.  They  found  that  the  boiling  heat  trans¬ 
fer  increases  and  decreases  for  Au-nanofluid  and  Si02-nanofluid, 
respectively.  These  conflicting  results  are  not  well-explained. 

Prakash  et  al.  [29,30]  performed  two  sets  of  experiments  to 
quantify  the  effect  of  heater  orientation  and  heater  surface  rough- 
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ness  on  pool  boiling  heat  transfer  of  A1203 /water-based  nanofluids. 
Their  first  set  of  experiment  [29]  with  vertical  tubular  heaters  of 
various  surface  roughness  showed  that  while  the  rough  heater 
surface  increases  heat  transfer,  smooth  surface  significantly  dete¬ 
riorates  the  heat  transfer.  For  example,  for  rough  heater  surface 
the  heat  transfer  enhancement  was  about  70%  at  0.5  wt.%  concen¬ 
tration  of  alumina  nanoparticles.  Whereas,  for  smooth  heater  the 
heat  flux  reduction  reaches  up  to  45%  at  a  particle  concentration 
of  2  wt.%.  Their  other  set  of  experiment  [30]  was  mainly  to  iden¬ 
tify  the  influence  of  heater  orientation  on  the  boiling  heat  transfer 
of  the  same  nanofluid  (Al203/water).  Their  results  demonstrated 
that  horizontal  and  inclined  heater  orientations  showed  enhance¬ 
ment  and  deterioration  of  boiling  heat  transfer  of  this  nanofluid, 
respectively. 

Soltani  et  al.  [31]  investigated  the  pool  boiling  heat  trans¬ 
fer  performance  of  Newtonian  nanofluids  under  various  heat  flux 
densities.  In  their  study,  7-Al203  (20-30  nm)/water  and  Sn02 
(55  nm)/water-based  nanofluids  were  used  in  a  vertical  cylin¬ 
drical  glass  vessel.  Their  results  showed  that  except  for  low 
concentrations  (>0.5  wt.%)  of  Sn02  nanoparticles,  the  boiling  heat 
transfer  coefficients  of  these  nanofluids  increase  with  increas¬ 
ing  concentration  of  nanoparticles.  These  paradoxical  results  were 
attributed  to  the  differences  in  thermal  conductivity  and  size  of 
these  two  nanoparticles.  In  a  similar  pool  boiling  experiment,  this 
same  research  group  [32]  recently  used  non-Newtonian  nanoflu¬ 
ids  of  7-Al203  nanoparticles  dispersed  in  carboxy  methyl  cellulose 
(cmc)-water  solution  and  the  boiling  heat  transfer  coefficient  was 
found  to  increase  by  about  25%  at  1 .4  wt.%  loading  of  nanoparticles. 
Recently,  Truong  et  al.  [33]  conducted  pool  boiling  experiments  of 
diamond,  ZnO,  and  A1203  nanoparticles-laden  aqueous  nanofluids 
with  modification  of  sandblasted  as  well  as  bare  plate  heaters.  They 
found  up  to  35%  increase  in  CHF  for  pre-coated  heaters  compared 
to  those  of  bare  and  sandblasted  heaters. 

Among  very  few  studies  on  flow  boiling  of  nanofluids,  Kim 
et  al.  [34]  reported  about  50%  enhancement  in  flow  boiling  CHF  for 
A1203 /water  nanofluids  flowing  through  a  vertical  stainless  steel 
tube.  Very  recently,  Henderson  et  al.  [35]  studied  refrigerant-based 
Si02  and  CuO-nanofluids  in  flow  boiling  experiments  in  horizon¬ 
tal  copper  tube.  They  found  that  while  the  boiling  heat  transfer 
coefficient  (BHTC)  of  Si02/R-134a  nanofluid  decreases  up  to  55%  in 
comparison  to  pure  R-134a,  the  BHTC  increases  more  than  100% 
for  CuO-laden  nanofluid  over  base  fluid,  i.e.  mixture  of  R-134a  and 
polyolester  oil  (PO). 

A  summary  of  studies  on  boiling  heat  transfer  of  nanofluids  is 
presented  in  Table  1.  It  can  be  noticed  that  besides  inconsistent 
and  contradictory  results,  most  of  the  researchers  used  alumina- 
nanofluids.  Whereas  very  few  efforts  have  been  made  with  carbon 
nanotubes  (CNTs)-nanofluids  which  exhibit  much  higher  thermal 
performance  compared  to  those  of  ceramic-nanofluids  [14,15].  It  is, 
therefore  worthwhile  to  study  boiling  heat  transfer  characteristics 
of  CNTs-nanofluids. 

A  comparison  of  heat  flux  versus  superheat  results  from  var¬ 
ious  groups  is  shown  in  Fig.  4.  From  these  representative  results 
(Fig.  4),  it  can  clearly  be  seen  that  heat  flux  (also  critical  heat  flux) 
data  relative  to  superheat  reported  by  various  research  groups  vary 
widely.  This  is  probably  due  to  the  differences  in  characterization  of 
nanofluids,  different  size  and  concentration  of  nanoparticles  used 
and  different  types  of  heaters  used  in  various  research  groups. 
Despite  of  the  fact  that  some  researchers  reported  deterioration 
of  boiling  heat  transfer  of  nanofluids,  the  significant  increase  in  the 
critical  heat  flux  of  nanofluids  is  still  undisputed. 

2.3.  Boiling  heat  transfer— correlations 

Although  there  are  many  correlations  [18,20]  to  predict 
the  boiling  heat  transfer  coefficient  and  critical  heat  flux, 
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Fig.  4.  Comparison  of  heat  flux  versus  superheat  results  from  various  research 
groups. 


researchers  working  on  boiling  heat  transfer  with  nanofluids 
[27,30,36,37]  mainly  used  predictions  from  two  most  popular 
correlations— Rohsenow  [45]  and  Zuber  [46]  correlations  to  com¬ 
pare  their  experimental  data.  Thus,  without  elaborating  details  of 
other  models  or  correlations,  these  two  classical  correlations  are 
presented  here. 

Based  on  dimensional  analysis  of  all  relevant  factors  arising  in 
nucleate  boiling  and  experimental  data  over  a  wide  variety  of  con¬ 
ditions  Rohsenow  [45]  proposed  one  of  the  first  and  most  widely 
used  nucleate  boiling  correlations  which  is  commonly  expressed  in 
the  following  form: 
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where  cp,  hfg,  Tw,  Ts,  q",  /x,  p,  pg,  a,gand  k  are  specific  heat  of  liquid, 
the  latent  heat  of  the  fluid,  temperature  of  heater  wall,  saturation 
temperature  of  liquid,  heat  flux,  viscosity  of  saturated  liquid,  den¬ 
sity  of  saturated  liquid,  density  of  saturated  vapor,  surface  tension, 
acceleration  of  gravity  and  thermal  conductivity  of  liquid,  respec¬ 
tively.  The  values  for  empirical  coefficient  Csj  and  exponent  n  for 
various  surface-fluid  combinations  can  be  found  elsewhere  [47]. 
For  example,  CSf  and  n  for  water-platinum  combination  are  0.013 
and  1,  respectively. 

The  maximum  heat  flux  attainable  in  the  nucleate  pool  boil¬ 
ing  regimes  is  the  critical  heat  flux  (CHF)  or  point  of  departure 
from  nucleate  boiling  (DNB).  A  widely  used  correlation  which  pre¬ 
dicts  the  pool  boiling  CHF  for  a  flat  infinite  heater  surface  (facing 
upwards)  is  Zuber’s  correlation  [46].  The  Zuber  CHF  is  expressed 
as: 


Qchf  =  KPg/2flfg  S^(P  -  pg) 


1  1/4 


The  value  of  constant  K  ranged  from  0.138  to  0.157.  However, 
simplifying  the  formulation  analysis,  Zuber  proposed  a  value  of 
I<=  0.131.  Although  widely  used,  most  of  the  cases  predictions  of 
pool  boiling  heat  transfer  and  critical  heat  flux  of  nanofluids  using 
these  two  abovementioned  correlations  do  not  agree  well  with  the 
experimental  results  [27,36].  There  is  therefore  a  strong  need  to 
develop  new  correlations  or  modify  these  classical  models  to  be 
used  for  nanofluids. 


2.4.  Present  boiling  experiments  with  CNT-nanofluids 

For  pool  boiling  experiments,  high  purity  single-walled  carbon 
nanotubes  (SWCNTs)  are  suspended  in  deionized  water  to  prepare 
sample  nanofluids.  As  a  part  of  surface  treatment,  CNT  bundles  were 
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Table  1 

Summary  of  boiling  experiments  with  nanofluids. 


Researchers 

Heater  (type  of  boiling) 

Nanofluids 

Remarks 

Yang  and  Maa  [7] 

Horizontal  tube  heater 
(pool) 

AI2O3 /water 

Boiling  heat  flux  increases  considerably. 

Witharana  [28] 

Cylindrical  vessels  (pool) 

Au  and  SiC>2 /water  and  EG 

Boiling  heat  transfer  increases  for  Au-nanofluids  but 
decreases  for  Si02-nanofluids. 

Das  et  al.  [26] 

Cylindrical  cartridge  (pool) 

AI2O3 /water 

Deterioration  of  heat  transfer  performance. 

You  et  al.  [24] 

Cartridge  (pool) 

AI2O3 /water 

Critical  heat  flux  (CHF)  increases  up  to  200%. 

Vassallo  et  al.  [36] 

NiCr  wire  (pool) 

SiC>2 /water 

CHF  increases  significantly. 

Bang  and  Chang  [37] 

Square  flat  heater  (pool) 

AI2O3 /water 

Pool  boiling  heat  transfer  deteriorates  but  CHF  increases. 

Wen  and  Ding  [27] 

Flat  disk  heater  (pool) 

AI2O3 /water 

The  boiling  heat  transfer  increases  up  to  40%  at  1 .2  wt.%  of 
nanoparticles. 

Kim  et  al.  [38] 

Smooth  NiCr  wire  heater 
(pool) 

TiC>2 /water 

The  maximum  enhancement  of  CHF  is  200%. 

Kim  et  al.  [39] 

Stainless  steel  wire  and  flat 
plate  heaters  (pool) 

AI2O3,  ZrC>2  and  SiC>2 /water 

The  maximum  increase  in  CHF  is  80%  but  the  heat  transfer 
rate  is  deteriorated. 

Jackson  [40] 

Flat  copper  coupon  heater 
(pool) 

Au/water 

While  the  heat  transfer  decreased  about  20%,  the 
maximum  CHF  increase  is  5  times  over  water. 

Prakash  et  al.  [29] 

Vertical  tubular  heaters 
(pool) 

AI2O3 /water 

While  rough  heater  surface  increases  heat  transfer,  smooth 
surface  significantly  deteriorates. 

Prakash  et  al.  [30] 

Tubular  heaters  at  various 
orientations  (pool) 

AI2O3 /water 

Horizontal  and  inclined  heater  orientations  showed 
enhancement  and  deterioration  of  heat  transfer, 
respectively. 

Chopkar  et  al.  [41  ] 

Flat  surface  (pool) 

ZrC>2 /water 

Enhanced  boiling  heat  transfer  is  found  at  low  particle 
concentration.  Tetramethyl  ammonium  hydroxide 
surfactant  increases  the  heat  transfer  rate. 

Liu  and  Liao  [42] 

Flat  horizontal  copper 
surface  (pool) 

CuO/water  and  Alcohol 
(C2H5OH) 

The  boiling  heat  transfer  characteristics  of  the  both 
nanofluids  are  somewhat  poor  and  the  CHF  values  are 
higher  than  those  of  the  base  liquids. 

Lv  and  Liu  [43] 

Vertical  small  heated  tubes 
(pool) 

CuO/water 

Enhancement  in  CHF  and  boiling  heat  transfer  rate  are 
observed  for  surfactant-free  nanofluids  but  addition 

NaDBS  surfactant  deteriorates  the  heat  transfer. 

Kathiravan  et  al.  [44] 

Horizontal  tube  (pool) 

CNTs/water 

Boiling  heat  transfer  coefficient  (BHTC)  increases  up  to 
1.75-fold  for  0.25  vol.%  of  CNT.  Sodium  lauryl  sulfate  (SLS) 
surfactant  of  9  wt.%  was  also  added. 

Kim  et  al.  [34] 

Vertical  stainless  steel  tube 
(flow) 

AI2O3 /water 

Flow  boiling  CHF  increases  about  50%. 

Soltani  et  al.  [31  ] 

Vertical  cylindrical  glass 
vessel  (pool) 

AI2O3  and  Sn02 /water 

Except  for  low  concentrations  (>0.5  wt.%)  of  SnC>2,  the 

BHTC  increases  with  loading  of  nanoparticles. 

Soltani  et  al.  [32] 

Vertical  cylindrical  glass 
vessel  (pool) 

AhOs/cmc-water  solution 

The  BHTC  increases  by  25%  for  1.4  wt.%  of  AI2O3 
nanoparticles. 

Henderson  et  al.  [35] 

Horizontal  copper  tube 
(flow) 

Si02/R-134a  and 

CuO/Rl  34a  +  polyolester  oil 

While  BHTC  of  SiC>2/R-134a  nanofluid  decreases  up  to  55% 
in  comparison  to  pure  R-134a,  it  increases  more  than  100% 
for  CuO-laden  nanofluid  over  base  fluid. 

Truong  et  al.  [33] 

Sandblasted  and  bare 
horizontal  plate  heaters 
(pool) 

Diamond,  ZnO  and 

AI2O3 /water 

The  CHF  of  the  pre-coated  heaters  increased  by  up  to  35% 
with  respect  to  that  of  the  bare,  sandblasted  heaters. 

refluxed  with  hydrochloric  acid  at  1 00  °C  for  several  hours.  This  acid 
was  chosen  as  a  reactive  reagent  because  it  removes  catalytic  par¬ 
ticles  without  reducing  the  length  of  the  tubes  or  damaging  the 
side  walls.  Other  acids  (HN03  and  HN03:H2S04)  tried  as  refluxing 
mixtures,  were  too  harsh  and  damaged  the  structure.  CNTs  were 
separated  from  the  acid  by  centrifugation  for  about  an  hour  and 
left  to  dry.  Different  concentrations  of  sodium  dodecyl  benzene- 
sulfonate  (NaDBS)  surfactant  are  used  as  dispersing  and  stabilizing 
agent  for  nanotubes  in  water.  The  schematic  of  experimental  facili¬ 
ties  used  is  shown  in  Fig.  5.  Details  of  the  experimental  facilities  and 
procedures  are  reported  in  previous  studies  [48,49]  and  will  not  be 
discussed  here.  Representative  results  of  surfactant-  and  surface 
tension-dependent  pool  boiling  experiments  of  CNTs-nanofluids 
[48]  are  presented  and  discussed  in  this  section.  Both  the  critical 
heat  flux  (CHF),  which  is  the  sudden  jump  in  temperature  at  one 
and  the  same  heat  flux  and  the  burnout  heat  flux  (BHF),  which  is  the 
case  when  the  wire  is  brought  to  catastrophic  failure  (i.e.  complete 
burnt  out)  are  determined  at  constant  0.1%  volumetric  loading  of 
CNT  and  for  various  concentrations  of  surfactant. 

2.4.2.  Effect  of  surfactant  concentration  on  boiling  heat  flux 

The  results  of  pool  boiling  heat  flux  with  respect  to  superheat 
(wire  temperature  minus  saturation  temperature  of  liquid,  Tw  -Ts) 
are  presented  and  compared  with  the  one  for  pure  de-ionized  water 


in  Fig.  6.  The  NaDBS  surfactant  to  CNTs  concentrations  are  varied 
from  1 :20  to  1 :1.  Fig.  6  demonstrates  that  the  CNT-nanofluid  with 
any  concentration  of  surfactant  exhibits  higher  CHF  value  than  that 
of  base  fluid.  The  CHF  value  of  pure  deionized  water  is  750  kW/m2. 
The  effect  of  increasing  the  surfactant  concentration  from  1 :20  to 
1 :5  (NaDBS:CNT)  is  that  the  CHF  increases  as  well.  However,  if  the 
concentration  of  surfactant  is  further  increased  from  1 :5  to  1 :1,  the 
CHF  drops  drastically  from  4439  kW/m2  to  1322  kW/m2,  respec¬ 
tively.  There  is,  therefore  a  critical  concentration  of  surfactant  for 
which  the  CHF  and  BHF  reach  maximum  values.  The  highest  CHF 
and  BHF  results  are  obtained  for  the  concentration  ratio  of  1:5. 
The  deposition  of  nanoparticles  on  the  heater  wire  is  believed  to 
be  one  of  the  key  reasons  for  any  enhancement  in  CHF  of  this 
nanofluid  up  to  critical  concentration  of  surfactant.  Kim  et  al.  [50] 
also  claimed  that  the  deposition  of  nanoparticles  on  the  heater  wire 
is  the  main  reason  for  the  enhancement  of  CHF  of  their  nanofluids. 
SEM  micrographs  shown  in  Fig.  7  demonstrate  that  with  an  increase 
in  surfactant  the  deposition  on  the  wire  increases  as  well.  However, 
the  maximum  CHF  and  BHF  do  not  occur  at  the  maximum  particle 
deposition,  i.e.  at  1:1  NaDBS:CNT  ratio. 

2.4.2.  Effect  of  surface  tension  on  burnout  heat  flux 

The  burnout  heat  flux  is  determined  at  constant  0.1%  volumetric 
loading  of  CNT  and  at  various  surface  tensions  of  nanofluids  which 
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Fig.  5.  Schematic  of  pool  boiling  experimental  setup  [48]. 
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Fig.  6.  Effect  of  surfactant  concentration  on  heat  flux  of  SWCNTs-nanofluid. 


was  adjusted  by  varying  concentrations  of  NaDBS  surfactant.  The 
surface  tensions  of  nanofluid  with  surfactant  (crNS)  and  deionized 
water  with  surfactant  (crws)  concentrations  are  measured.  At  any 
concentration  of  surfactant,  the  difference  between  surface  ten¬ 
sions  of  nanofluid  and  its  based  fluid  (deionized  water)  is  termed 
as  relaxation  of  surface  tension,  i.e.  Aar  =  crNS  -  orWs-  This  relaxation 
of  surface  tension  of  nanofluids  is  the  driving  force  for  the  burnout 
heat  flux  of  nanofluid,  i.e.  (BHF)N  in  pool  boiling  as  shown  in  Fig.  8. 
Based  on  the  best  fit  of  the  experimental  data,  a  nonlinear  empir¬ 
ical  correlation  between  the  BHF  and  relaxation  of  surface  tension 


G  -a  (mN/m) 

NS  WS 


Fig.  8.  Effect  of  surface  tension  on  the  burnout  heat  flux  of  nanofluid  (NaDBS  to 
CNTs  concentration  ratios  are  shown  in  parentheses). 


is  obtained  and  has  the  form 

(BHF)N  =  1  +  °'073 A<7r  -  0.0065Aar2  +  0.0003 Aar3 


where  (BF1F)W  is  the  burnout  heat  flux  of  water  and 

Acrr  =  (ctns  -  ^ws)  mN/m. 

It  is  demonstrated  that  regardless  of  surfactant  concentration, 
the  burnout  heat  flux  of  nanofluid  strongly  depends  (i.e.  nonlinearly 
related)  on  its  relaxation  of  surface  tension  with  base  fluid.  Similar 
to  the  CHF,  the  BF1F  of  nanofluid  also  reaches  its  highest  value  at 
NaDBS  and  CNT  concentration  ratio  of  1 :5.  The  highest  BHF  value  of 
nanofluid  is  about  265%  enhancement  over  that  for  the  base  fluid  at 
A  or «  25  mN/m.  The  burnout  heat  flux  of  deionized  water  (BHF)W  is 
found  to  be  1 500  kW/m2 .  The  observed  increment  of  relative  BHF  of 
nanofluid  (i.e.  (BHF)N/(BHF)W)  can  be  attributed  to  the  Marangoni 
effect  which  results  from  the  surface  tension  differences.  Both  the 
CHF  and  burnout  are  found  to  happen  at  the  optimum  1:5  sur- 
factant:CNT  ratio.  Thus,  the  pool  boiling  heat  transfer  behavior  of 
surfactant-added  nanofluid  is  not  only  dictated  by  the  deposition 
of  nanotubes  on  heater  surface,  but  also  by  the  relaxation  of  sur¬ 
face  tension  of  nanofluid,  which  is  a  precursor  to  the  deposition  and 
overall  onsets  of  both  the  CHF  and  BHF. 


2.5.  Studies  on  nanofluids  droplet  spreading 

Droplets  impingement  and  spreading  studies  on  heated  sub¬ 
strate  surface  are  of  great  importance  for  the  practical  application 
of  nanofluids  as  an  advanced  coolant  particularly  in  spray  cooling. 
Unfortunately,  very  little  research  efforts  have  been  made  to  study 
nanofluids’  droplet  impingement  behavior  on  solid  surfaces  under 
various  conditions.  Among  very  few  studies,  Wasan  and  Nikolov 
[51  ]  were  the  first  to  investigate  the  effects  of  the  particle  structure 
formation  and  the  structural  disjoining  pressure  of  nanoparti¬ 
cles  on  the  spreading  of  nanofluids  on  solid  surface.  Recently, 
Duursma  et  al.  [52]  studied  the  effect  of  aluminum  nanoparti¬ 
cles  on  droplets  boil-off  by  allowing  nanofluid  drops  to  fall  onto 
a  copper  surface  at  temperature  higher  than  the  liquid  satura- 
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Fig.  7.  SEM  images  of  NiCr  wires  at  various  surfactant  concentrations  after  pool  boiling  experiments  [48]. 
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tion  temperature.  They  demonstrated  that  increasing  the  surface 
temperature  and  Weber  number  promote  the  receding  breakup 
scenario,  while  increasing  the  nanoparticle  concentration  discour¬ 
ages  this  breakup.  In  another  recent  study,  the  influence  of  surface 
temperature  on  the  hydrodynamic  characteristics  of  water  and  a 
nanofluid  droplets  impinging  on  a  polished  and  nanostructured 
surface  was  investigated  by  Shen  et  al.  [53].  Their  results  showed 
that  SWCNTs-nanofluid  has  larger  spreading  diameter  compared  to 
that  of  deionized  water  and  using  a  nanofluid  or  a  nanostructured 
surface  can  reduce  the  total  evaporation  time  up  to  37%.  Murshed 
and  Nieto  de  Castro  [54]  very  recently  reported  the  spreading  char¬ 
acteristics  of  a  nanofluid  droplets  impinging  on  a  metallic  substrate 
under  the  influence  of  several  key  factors  such  as  nanoparticles 
volume  fraction,  substrate  temperature,  and  the  Weber  number. 
Results  showed  that  the  transient  spreading  diameter  and  height 
of  droplet  impacting  onto  solid  surface  are  greatly  influenced  by 
each  of  these  factors.  Nevertheless,  more  studies  are  needed  on 
dynamics  of  non-boiling  and  boiling  nanofluid  droplets  impinging 
on  solid  or  liquid  surfaces  as  the  spreading  of  liquid  droplet  plays 
a  key  role  in  many  industrial  processes  like  spray  cooling,  coating, 
ink-jet  printing,  and  oily  soil  removal. 


3.  Flow  and  heat  transfer  characteristics  of  nanofluids 

Compared  to  the  reported  studies  on  thermal  conductiv¬ 
ity,  investigations  on  convective  heat  transfer  of  nanofluids  are 
still  scarce.  However,  the  practical  applications  of  nanofluids  as 
advanced  heat  transfer  fluids  are  mainly  in  flowing  systems  such  as 
mini-  or  micro-channels  heat  sinks  and  miniaturized  heat  exchang¬ 
ers.  In  this  section,  we  critically  review  both  experimental  and 
computational  studies  on  convective  heat  transfer  of  nanofluids. 
In  addition,  some  representative  results  from  our  experimental 
investigation  on  laminar  flow  convective  heat  transfer  of  Ti02- 
nanofluids  are  also  presented  and  discussed. 


3.1.  Experimental  studies  on  convective  heat  transfer  of 
nanofluids 

The  first  experiment  on  convective  heat  transfer  of  nanofluids 
(e.g.  7-AI2O3 /water)  under  turbulent  flow  conditions  was  per¬ 
formed  by  Pak  and  Cho  [55].  In  their  study,  even  though  the  Nusselt 
number  (Afti)  was  found  to  increase  with  increasing  nanoparticle 
volume  fraction  and  Reynolds  number,  the  heat  transfer  coefficient 
(ft)  actually  decreased  by  3-12%.  The  reasons  for  such  paradoxical 
results  might  be  the  observed  large  enhancement  in  viscosity.  On 
the  other  hand,  Eastman  et  al.  [56]  later  showed  that  with  less  than 
1  vol.%  of  CuO  nanoparticles,  the  convective  heat  transfer  coeffi¬ 
cient  (ft)  of  water  increased  more  than  15%.  The  results  of  Xuan 
and  Li  [57]  illustrated  that  the  Nusselt  number  of  Cu/water-based 
nanofluids  increased  significantly  with  the  volumetric  loading 
of  particles  as  for  2  vol.%  of  nanoparticles,  the  Nusselt  number 
increased  by  about  60%.  Wen  and  Ding  [58]  investigated  the  heat 
transfer  behavior  of  nanofluids  at  the  tube  entrance  region  under 
laminar  flow  conditions  and  showed  that  the  local  heat  transfer 
coefficient  varied  with  particle  volume  fraction  (0)  and  Reynolds 
number  {Re).  They  also  observed  that  the  enhancement  is  par¬ 
ticularly  significant  at  the  entrance  region.  Later  Heris  et  al.  [59] 
studied  convective  heat  transfer  of  CuO  and  A1203 /water-based 
nanofluids  under  laminar  flow  conditions  through  an  annular  tube. 
Their  results  showed  that  heat  transfer  coefficient  increases  with 
particle  volume  fraction  as  well  as  Peclet  number.  In  their  study, 
A1203 /water-based  nanofluids  found  to  have  higher  enhancement 
of  heat  transfer  coefficient  compared  to  CuO/water-based  nanoflu¬ 
ids. 
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Fig.  9.  Comparison  convective  heat  transfer  results  from  various  research  groups 
[55,57,58,62,70]. 

An  experimental  investigation  on  the  forced  convective  heat 
transfer  and  flow  characteristics  of  Ti02 -water  nanofluids  under 
turbulent  flow  conditions  is  reported  by  Duangthongsuk  and 
Wongwises  [60].  A  horizontal  double-tube  counter  flow  heat 
exchanger  is  used  in  their  study.  They  observed  a  slightly  higher 
(6-11%)  heat  transfer  coefficient  for  nanofluid  compared  to  pure 
water.  The  heat  transfer  coefficient  increases  with  increasing  mass 
flow  rate  of  the  hot  water  as  well  as  nanofluid.  They  also  claimed 
that  the  use  of  the  nanofluid  has  a  little  penalty  in  pressure  drop. 

In  microchannel  flow  of  nanofluids,  Faulkner  et  al.  [61  ]  was  the 
first  to  conduct  convective  heat  transfer  experiments  with  aque¬ 
ous  CNT-nanofluid  in  a  microchannel  with  hydraulic  diameter  of 
355  |jim  at  Reynolds  numbers  between  2  and  17.  They  found  con¬ 
siderable  enhancement  in  heat  transfer  coefficient  of  this  nanofluid 
at  CNT  concentration  of  4.4%.  Later,  Jung  et  al.  [62]  studied  heat 
transfer  performance  of  A1203 /water-based  nanofluid  in  a  rectan¬ 
gular  microchannel  under  laminar  flow  condition.  Results  showed 
that  the  heat  transfer  coefficient  increased  by  more  than  32%  for 
1.8  vol.%  of  nanoparticles  and  the  Nusselt  number  increases  with 
increasing  Reynolds  number  in  the  flow  regime  of  5  >  Re  <  300. 

An  up-to-date  overview  of  the  published  experimental  results 
on  the  convective  heat  transfer  characteristics  of  nanofluids  is 
given  in  Table  2.  A  comparison  of  results  of  Nusselt  number  versus 
Reynolds  number  for  both  laminar  and  turbulent  flow  conditions 
from  various  groups  is  shown  in  Fig.  9.  From  Table  2  and  Fig.  9,  it 
can  be  seen  that  the  results  from  various  groups  vary  widely  and 
most  of  the  studies  lack  physics-based  explanation  of  the  observed 
results.  Although  there  are  couple  of  review  articles  [71,72]  that 
attempted  to  put  together  some  studies  on  convective  heat  trans¬ 
fer  of  nanofluids,  no  comprehensive  analysis  of  up-to-date  findings 
are  reported. 

Among  very  few  research  efforts,  the  work  of  Putra  et  al.  [73] 
was  the  first  to  investigate  the  natural  convective  heat  transfer  of 
aqueous  CuO  and  Al203-nanofluids.  They  used  a  horizontal  poly- 
oxymethylene  cylinder  which  was  heated  from  one  end  and  cooled 
from  the  other.  Significant  deterioration  (decrease)  of  convection 
heat  transfer  for  these  nanofluids  was  observed  and  the  deteriora¬ 
tion  increases  with  increasing  particle  concentration,  particularly 
for  CuO  nanoparticles.  Putra  et  al.  ascribed  the  possible  reasons  of 
such  deterioration  to  the  effects  of  particle-fluid  slip  and  sedimen¬ 
tation  of  nanoparticles.  In  contrast  to  Putra  et  al.  [73],  numerical 
simulation  of  Khanafer  et  al.  [74]  showed  that  in  a  2D  horizon¬ 
tal  enclosure,  the  natural  convection  heat  transfer  coefficient  (Afti) 
of  nanofluids  increases  with  particle  concentration.  By  taking  into 
account  the  solid  particle  dispersion,  they  [74]  developed  a  model 
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Table  2 

Summary  of  forced  convection  heat  transfer  experimental  studies  of  nanofluids. 
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Researchers 

Geometry/flow  nature 

Nanofluids 

Findings 

Pak  and  Cho  [55] 

Tube/turbulent 

AI2O3  and  TiC>2 /water 

At  3  vol  %,  h  was  12%  smaller  than  pure  water  for  a  given 
average  fluid  velocity. 

Xuan  and  Li  [57] 

Tube/turbulent 

Cu/water 

A  larger  enhancement  of  h  with  volume  fraction  (0)  and  Re 
was  observed. 

Wen  and  Ding  [58] 

Tube/laminar 

AI2O3 /water 

Increased  h  with  0  and  Reynolds  number  was  observed. 
NaDBS  was  used  as  the  surfactant 

Ding  et  al.  [63] 

Tube/laminar 

CNT/water 

At  0.5  wt.%,  h  increased  by  more  than  350%  at  Re  =  800. 

Yang  et  al.  [64] 

Tube/laminar 

Graphite/automatic 
transmission  fluid 

The  nanoparticles  considerably  increase  the  heat  transfer 
coefficient  of  the  fluid  system  in  laminar  flow. 

Heris  et  al.  [59] 

Tube/laminar 

AI2O3  and  CuO/water 

h  increase  with  0  and  Pe.  AI2O3  shows  higher 
enhancement  than  that  of  CuO. 

Lai  et  al.  [65] 

Tube/laminar 

AI2O3 /water 

Nu  increased  8%  for  0  =  0.01  and  Re  =  270. 

Jung  et  al.  [62] 

Rectangular 

microchannel/laminar 

A1203 /water 

h  increased  1 5%  for  0  =  0.01 8. 

Williams  et  al.  [66] 

Tube/turbulent 

AI2O3  and  ZrC>2 /water 

h  increased  significantly. 

Hwang  et  al.  [67] 

Tube/laminar 

A1203 /water 

h  increased  only  up  to  8%  at  Re  =  730  for  0  =  0.003. 

Xie  et  al.  [68] 

Tube/laminar 

AI2O3,  ZnO,  Ti02  and 
MgO/water 

h  increased  up  to  252%  at  Re  =  1000  for  MgO  nanofluid. 

Amrollahi  et  al.  [69] 

Tube/laminar  and 
turbulent 

MWCNT/water 

h  increased  up  to  33-40%  at  concentration  of  0.25  wt.%. 

to  analyze  the  heat  transfer  performance  of  nanofluids  inside  an 
enclosure. 

Wen  and  Ding  [75]  later  investigated  heat  transfer  behavior  of 
specially  formulated  Ti02 /water  nanofluid  under  the  natural  con¬ 
vection  conditions.  The  results  showed  that  the  natural  convective 
heat  transfer  coefficient  decreases  with  increasing  particle  con¬ 
centration.  These  unexpected  results  are  in  contradiction  to  the 
numerical  findings  of  Khanafer  et  al.  [74]  but  are  in  agreement  with 
the  observations  by  Putra  et  al.  [73  ].  Along  with  the  reasons  ascribed 
by  Putra  et  al.  [73],  Wen  and  Ding  [75]  added  several  other  rea¬ 
sons  for  such  paradoxical  results  that  include  convection  induced 
by  concentration  difference,  modification  of  dispersion  properties 
and  particle-surface  and  particle-particle  interactions. 

3.2.  Theoretical  models  for  convective  heat  transfer  of  nanofluids 

Researchers  investigating  convective  heat  transfer  of  nanoflu¬ 
ids  employed  existing  conventional  single-phase  fluid  correlation 
to  predict  the  heat  transfer  coefficient  or  proposed  new  correlations 
obtained  by  fitting  their  own  experimental  data.  However,  none  of 
these  correlations  were  validated  with  wide  range  of  experimental 
data  under  various  conditions  and  thus  are  not  widely  accepted.  In 
an  attempt  to  establish  a  strong  explanation  of  the  reported  anoma¬ 
lously  enhanced  convective  heat  transfer  coefficient  of  nanofluid, 
Buongiorno  [76]  considered  seven-slip  mechanisms  and  concluded 
that  among  those  seven  only  Brownian  diffusion  and  thermophore¬ 
sis  are  the  two  most  important  particle/fluid  slip  mechanisms  in 
nanofluids.  Besides  proposing  a  new  correlation,  he  also  claimed 
that  the  enhanced  laminar  flow  convective  heat  transfer  can  be 
attributed  to  a  reduction  of  viscosity  within  and  consequent  thin¬ 
ning  of  the  laminar  sublayer.  Commonly  used  classical  models 
together  with  recently  proposed  correlations  for  the  predictions 
of  convective  heat  transfer  coefficient  of  nanofluids  are  summa¬ 
rized  in  Table  3.  Considering  all  possible  key  mechanisms  such  as 
particle  Brownian  motion  and  migration,  thermophoresis  and  iner¬ 
tial  a  generalized  convective  heat  transfer  correlation  needs  to  be 
developed. 

3.3.  Numerical  studies  on  convective  heat  transfer  of  nanofluids 

Several  preliminary  numerical  studies  on  the  convective  heat 
transport  of  nanofluids  were  conducted  by  Khanafer  et  al.  [74] 
and  a  research  group  from  Universite  de  Moncton,  Canada  [79,80]. 
Results  from  the  Canadian  group  showed  that  nanofluids  can  sub¬ 


stantially  enhance  heat  removal  in  both  tube  [79]  and  radial  flows 
[80].  For  example,  Roy  et  al.  [80]  numerically  studied  the  hydro- 
dynamic  and  thermal  fields  of  A1203 /water-based  nanofluids  in  a 
radial  laminar  flow  system  and  reported  a  two-fold  increase  in  heat 
transfer  for  10vol.%  of  nanoparticles  in  water.  They  also  observed 
that  the  wall  shear  stress  increased  with  an  increase  in  particle  vol¬ 
ume  fraction.  For  Al203/water-  and  Al203/ethylene  glycol-based 
nanofluids,  numerical  results  of  Mai'ga  et  al.  [79]  showed  that  the 
heat  transfer  coefficient  increases  considerably  with  increasing 
nanoparticle  concentration.  However,  the  presence  of  such  parti¬ 
cles  has  also  induced  drastic  effects  on  the  wall  shear  stress  that 
increases  appreciably  with  the  particle  loadings. 

A  Lattice  Boltzmann  model  for  simulating  flow  and  energy  trans¬ 
port  processes  inside  nanofluids  was  proposed  by  Xuan  et  al.  [81  ]. 
They  considered  the  external  and  internal  forces  acting  on  the  sus¬ 
pended  nanoparticles  and  interactions  among  the  nanoparticles 
and  fluid  particles.  Based  on  the  simulation  results,  they  inferred 
that  the  fluctuation  of  Nusselt  number  of  nanofluids  along  the  main 
flow  direction  is  due  to  the  unstable  distribution  of  nanoparticles. 
No  further  development  is  made  in  exploring  the  use  of  the  Lattice 
Boltzmann  method  to  explain  the  anomalous  thermal  convection 
of  nanofluids. 

Palm  et  al.  [82]  performed  numerical  simulations  for  the  laminar 
forced  convection  flow  of  nanofluids  with  temperature-dependent 
properties.  They  found  that  A1203 /water-based  nanofluids  with 
a  volume  fraction  of  4%  can  produce  a  25%  increase  in  the 
average  wall  heat  transfer  coefficient  compared  to  water  alone. 
Significant  differences  were  also  found  when  using  constant 
property  nanofluids  (temperature-independent)  versus  nanofluids 
with  temperature-dependent  properties.  The  use  of  temperature- 
dependent  properties  resulted  in  larger  heat  transfer  predictions 
with  corresponding  decrease  in  wall  shear  stresses  when  com¬ 
pared  to  predictions  using  constant  properties  nanofluids.  With  an 
increase  in  wall  heat  flux,  they  found  that  the  average  heat  transfer 
coefficient  increases  while  the  wall  shear  stress  decreases. 

Based  on  Khanafer  et  al.’s  [74]  model  and  taking  into  account 
particle  dispersion,  Jou  and  Tzeng  [83]  numerically  studied 
convective  heat  transfer  performance  of  nanofluids  inside  two- 
dimensional  rectangular  enclosures.  Their  results  showed  that 
increasing  volume  fraction  and  buoyancy  parameter  can  cause  an 
increase  in  average  heat  transfer  coefficient. 

The  cooling  performance  of  a  microchannel  (silicon)  heat  sink 
under  forced  convective  laminar  flow  with  nanofluids  is  numeri¬ 
cally  investigated  by  Jang  and  Choi  [84].  Results  showed  that  at  a 
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Table  3 

Correlations  used  for  predicting  the  convective  heat  transfer  coefficient  of  nanofluids. 


Origin/reference 


Correlations 


Flow  regimes/conditions/remarks 


Shah  equation  [77] 


1.953  (  Re  Pr- 
x 


1/3 


D 


for  (  Re Pr^  ]  >  33.3 
D 


4.364  +  0.0722 RePr-;  for  Re Pr-  <  33.3 


x 


x 


Laminar  flow  and  constant  heat  flux  condition.  It  is  popular 
for  a  thermal  entrance  region. 


Dittus  and  Boelter  [78] 

Nu  =  0.023 Re0-8  Pr0-4 

Heating  of  fluid  and  used  for  0.7  <Pr<  100.  It  is  very 
popular  in  pipe  flow  cases  and  is  not  recommended  for 
temperature  dependent  property  variation. 

Sieder-Tate  equation  [77] 

Nu  =  0.027 Re0-8  Pr^ 13  (iib  1 /jls)0A4 

Valid  for  0.7  <Pr<  16,700  and  Re>  10,000.  It  is  used  for 

where  /ib  and  /xs  are  the  fluid  viscosity  at  bulk 
fluid  temperature  and  at  heat-transfer  boundary 
surface  temperature,  respectively. 

temperature-dependent  property  variation  cases. 

Pak  and  Cho  [55] 

Nu  =  0.021  Re0-8  Pr0-5 

Turbulent  flow 

Xuan  and  Li  [57] 

Laminar  flow: 

Pe  is  the  particle  Peclet  number.  It  considered  particles 

Nu  =  0.4328(1  + 1 1 .285(p0J54Pe0218)Re033Pr0A 

Turbulent  flow: 

Nu  =  0.0059(1  +  7.62860a6886Pe0001  )Pe0-9238  Pr0A 

concentration. 

Mai'ga  et  al.  [79] 

Averaged  Nu  for  constant  heat  flux: 

Correlations  are  proposed  from  numerical  study  and  valid 

Nu  =  0.086 Re°-55Pr0-5 

for  Pe  <  1 000  and  6  <  Pr  <  753. 

Averaged  Nu  for  constant  temperature: 

Nu  =  0.28  Re°-35Pr036 

Jung  et  al.  [62] 

Nu  =  0.0095 0001  Re0A2Pr°-55 

For  microchannel  laminar  flow  regime.  It  considered  the 
volume  fraction  of  nanoparticles. 

Buongiorno  [76] 

N  (f/8)(Reb-\000)Prb 

l+5+y?/8(P+/3-l) 

Turbulent  flow 

where/is  the  friction  factor  and  SA  is  an  empirical 

constant. 

fixed  pumping  power  (2.245  W)  the  cooling  performance  of  their 
heat  sink  with  aqueous-diamond  (1  vol.%)  nanofluid  is  enhanced 
by  about  10%  compared  to  that  of  with  water  alone.  They  also 
reported  that  nanofluids  reduced  both  the  thermal  resistance  and 
the  temperature  difference  between  the  heated  microchannel  wall 
and  coolant. 

Using  a  dispersion  model  Heris  et  al.  [85]  later  performed 
numerical  simulations  for  the  laminar-flow  convective  heat  trans¬ 
fer  of  aqueous  nanofluids  in  a  circular  tube  with  constant  wall 
temperature.  Effects  of  concentration  and  size  of  various  types 
of  nanoparticles  (A1203,  CuO,  and  Cu  in  water)  as  well  as  Peclet 
number  on  heat  transfer  coefficient  are  investigated.  They  reported 
that  Nu  increases  with  increasing  particle  concentration  and  with 
decreasing  nanoparticles  size.  Their  numerical  predictions  were 
found  to  be  in  considerable  agreement  with  their  experimental 
data. 

A  numerical  study  on  laminar  mixed  convection  of  Al203/water 
nanofluid  was  reported  by  Mirmasoumi  and  Behzadmehr  [86]. 
Two-phase  mixture  model  was  used  to  study  the  effects  of  nanopar¬ 
ticles  mean  diameter  on  the  flow  parameters.  Their  calculated 
results  showed  that  like  Heris  et  al.’s  [85]  findings,  the  heat  transfer 
coefficient  significantly  increases  with  decreasing  the  nanoparti¬ 
cles  mean  diameter.  However  particle  size  does  not  notably  change 
the  skin  friction  coefficient. 

In  a  separate  numerical  study,  Ho  et  al.  [87]  attempted  to  identify 
the  effects  of  uncertainties  in  various  effective  viscosity  and  ther¬ 
mal  conductivity  models  adopted  on  buoyancy-driven  convection 
heat  transfer  of  alumina-water  nanofluid  in  a  square  enclosure. 
Based  on  their  simulation  results  from  a  comparative  study  of  four 
different  models,  they  demonstrated  that  the  effect  of  nanopar¬ 
ticle  concentration  on  the  averaged  Nusselt  number  significantly 
depends  on  the  effective  viscosity  and  thermal  conductivity  models 
used. 

Oztop  and  Abu-Nada  [88]  recently  numerically  investigated 
the  buoyancy-driven  natural  convection  heat  transfer  character¬ 
istics  of  different  types  of  water-based  nanofluids  containing  Cu, 
A1203,  and  Ti02  nanoparticles  in  a  partially  heated  enclosure.  An 
increase  in  mean  Nusselt  number  was  observed  with  the  volume 


fraction  of  nanoparticles  for  the  whole  range  of  Rayleigh  number 
( 1 03  <  Ra  <  5  x  1 05 )  and  heat  transfer  also  increases  with  increasing 
of  height  of  heater. 

A  very  recent  numerical  work  on  a  3D  laminar  flow  heat  trans¬ 
fer  with  A1203,  CuO/mixture  of  ethylene  glycol  and  water-based 
nanofluids  circulating  through  flat  tubes  of  a  radiator  was  under¬ 
taken  by  Vajjha  et  al.  [89].  Their  results  demonstrated  that  the 
convective  heat  transfer  coefficient  in  the  developing  and  devel¬ 
oped  regions  along  the  flat  tubes  with  the  nanofluid  flow  showed 
remarkable  improvement  over  the  base  fluid  and  for  the  local 
and  the  average  friction  factor  and  convective  heat  transfer  coef¬ 
ficient  show  an  increase  with  increasing  volumetric  concentration 
of  nanoparticles. 

Above  review  on  numerical  studies  demonstrated  that  convec¬ 
tive  heat  transfer  coefficient  of  nanofluids  considerably  increases 
with  increasing  concentration  of  nanoparticles  and  with  decreas¬ 
ing  nanoparticle  size.  However,  more  works  on  forced  and  natural 
convective  heat  transfer  of  nanofluids  are  to  be  performed  by  using 
advanced  numerical  methods  such  as  Lattice  Boltzmann,  Molecular 
and  Brownian  dynamic  simulations. 

3.4.  Present  experiments  on  laminar  flow  convective  heat  transfer 
of  nanofluids 

Sample  nanofluids  for  convective  heat  transfer  experiments 
were  prepared  by  dispersing  different  volume  percentages,  i.e. 
0.2-0.8%  of  titanium  dioxide  (Ti02 )  nanoparticles  ( 1 5  nm  diameter) 
in  deionized  water  (DIW).  As  a  dispersant  agent,  Cetyl  Trimethyl 
Ammonium  Bromide  (CTAB)  surfactant  was  added  in  the  sample. 
To  ensure  proper  dispersion  of  nanoparticles,  sample  nanofluid  was 
homogenized  by  using  an  ultrasonic  dismembrator. 

The  effects  of  nanoparticles  concentration  and  Reynolds  number 
on  the  convective  heat  transfer  coefficient  of  Ti02/deionized  water 
based-nanofluids  were  previously  investigated  [90]  and  some  rep¬ 
resentative  results  are  presented.  An  experimental  setup  was 
established  to  conduct  experiments  on  heat  transfer  of  nanofluids 
at  laminar  flow  regime  in  a  cylindrical  channel.  The  schematic  of 
experimental  facilities  used  is  shown  in  Fig.  1 0.  Details  of  the  exper- 
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Fig.  10.  Schematic  of  convective  heat  transfer  experimental  setup  [90]. 


imental  facilities  and  procedures  are  reported  elsewhere  [90]  and 
will  not  be  discussed  here.  Instead  formulations  for  obtaining  the 
heat  transfer  coefficient  (ft)  together  with  a  most  commonly  used 
existing  correlation  for  predicting  the  Nusselt  numbers  of  nanoflu¬ 
ids  in  laminar  flow  conditions  are  presented. 

As  detailed  in  previous  paper  [90],  applying  first  law  (energy  bal¬ 
ance)  in  control  volume  of  flow  channel,  the  following  formulation 
for  the  local  heat  transfer  coefficient  is  obtained 

Kf-x  = 

_ of _ 

{To,wW  ~  (fa[2 D20  ln(D0/D,9  -  (D2  -  D2)])/4tt(D20  -  D2)ksx)}  -  [Tt  +  ((T0  -  T,)/L)x} 

(4) 


where  T0?w(x)  is  the  outer  wall  temperature  of  the  tube  (mea¬ 
surable),  q"  is  the  heat  flux  of  the  test  section  (W/m2),  q  is  the  heat 
supplied  to  the  test  section  (W),  ks  is  the  thermal  conductivity  of  the 
copper  tube  (W/mK),  D,  and  D0  are  the  inner  and  outer  diameters 
of  the  tube,  respectively  and  x  represents  the  longitudinal  location 
of  the  section  of  interest  from  the  entrance.  L  is  the  length  of  the 
test  section  and  f  and  T0  are  the  inlet  and  outlet  fluid  temperature, 
respectively. 

Once  the  local  heat  transfer  coefficient  is  determined  and  the 
thermal  conductivity  of  the  medium  is  known,  the  local  Nusselt 
number  is  calculated  from 


NUnf-x 


huf-x^i 

k nf 


where  knf  is  the  effective  thermal  conductivity  of  nanofluids.  The 
classical  Hamilton-Crosser  model  is  used  for  the  determination  of 
knf  which  is  given  by  [91  ] 


kp  +  (n  -  1  )kf  -  (n  -  1  )(p{kf  -  1<p ) 
kp  +  (n  —  1  )kj  +  0(kj  —  kp) 


where  kf  and  kp  are  the  thermal  conductivities  of  the  base  liquid  and 
the  nanoparticles,  respectively,  0  is  the  volume  fraction  of  nanopar¬ 
ticles  and  n  is  the  empirical  shape  factor,  which  has  a  value  of  3  for 
spherical  particle. 

The  Nusselt  number  can  also  be  determined  from  the  existing 
correlations.  The  well-known  Shah’s  correlation  for  laminar  flows 
under  the  constant  heat  flux  boundary  conditions  is  used  and  repro¬ 
duced  as [77] 


/  D-X1/3  /  D\ 

Nu  =  1.953  ( RePr-^  J  for  iRePr^j  >33.3 
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Fig.  11.  Axial  profiles  of  local  heat  transfer  coefficient  of  nanofluid  at  Re  =  1100  [90]. 


For  steady  and  incompressible  flow  of  nanofluids  in  a  tube  of 
uniform  cross-sectional  area,  the  Reynolds  number  and  Prandtl 
number  are  defined  as  follows 


Re  = 


4m 


TlDi/lnf 


and  Pr  = 


-p-n/Mn/ 

Kf 


where  m  is  the  mass  flow  rate  and  /inf  and  cp_nf  are  the  viscosity 
and  specific  heat  of  nanofluids,  respectively. 

While  the  specific  heat  of  nanofluids  is  calculated  using  the  fol¬ 
lowing  volume  fractioned-based  mixture  rule  [55,79,84] 


Cp-nf  —  tpCp-p  T  (1  (9) 

the  viscosity  of  nanofluids  is  determined  from  Batchelor’s  model 
given  by  [92] 

/inf  =  M/(l  +  2-50  +  6.2  02)  (10) 

where  0  is  particle  volume  fraction  and  /if  is  the  base  fluid  viscosity. 
It  is  noted  that  other  classical  models  for  calculating  the  viscosity 
of  mixture  also  yield  similar  results  [13]. 


3.4 A.  Axial  profiles  of  local  heat  transfer  coefficient 

Fig.  1 1  illustrates  the  local  heat  transfer  coefficient  against  the 
axial  distance  from  the  entrance  of  the  test  section  at  Reynolds 
number  (Re)  of  1 1 00.  The  results  show  that  nanofluids  exhibits  con¬ 
siderably  enhanced  convective  heat  transfer  coefficient  which  also 
increases  with  volumetric  loadings  of  Ti02  nanoparticles.  For  exam¬ 
ple,  at  0.8vol.%  of  nanoparticles  and  at  position  x/D/ =  25  (where 
tube  diameter  Di  =  4mm),  the  local  heat  transfer  coefficient  of  this 
nanofluid  was  found  to  be  about  12%  higher  compared  to  deionized 
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Fig.  12.  Comparison  with  Shah’s  correlation  along  axial  distance  at  Re  =  1100  and 
1700  [90]. 
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Fig.  13.  Effect  of  Reynolds  number  and  nanoparticle  volume  fraction  on  Nusselt 
number  at  axial  location  of  x/D,  =  25  [90]. 
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water  at  the  same  Re.  The  enhancement  in  heat  transfer  coefficients 
of  nanofluids  is  believed  due  to  the  enhanced  effective  thermal  con¬ 
ductivity  and  the  acceleration  of  the  energy  exchange  process  in  the 
fluid  due  to  the  random  movements  of  the  nanoparticles.  Another 
reason  for  such  enhancement  can  be  the  migration  of  nanoparticles 
in  base  fluids  due  to  shear  action,  viscosity  gradient  and  Brownian 
motion  in  the  cross  section  of  the  tube. 

Fig.  12  compares  experimentally  determined  Nusselt  numbers 
with  the  predictions  by  Shah’s  correlation,  i.e.  Eq.  (7)  along  the 
axial  distance.  Shah’s  correlation  slightly  over-predicts  the  present 
results.  The  difference  in  tube  size  may  be  one  of  the  reasons  for 
such  over  prediction.  A  relatively  small  tube  (4  mm  diameter)  was 
used  in  this  experiment,  whereas  the  Shah’s  equation  was  devel¬ 
oped  on  the  basis  of  laminar  flow  in  large  channels.  Nevertheless, 
similar  over  prediction  by  Shah’s  equation  was  also  reported  by 
Wen  and  Ding  [58]. 

3.4.2.  Effect  of  Reynolds  number  and  particle  volume  fraction  on 
Nusselt  number 

The  effect  of  Reynolds  number  on  Nusselt  number  is  shown 
in  Fig.  13.  It  can  be  seen  that  the  measured  Nusselt  numbers  for 
nanofluids  are  higher  than  those  of  water  and  they  increase  remark¬ 
ably  with  Reynolds  number.  The  observed  enhancement  of  the 
Nusselt  number  could  be  due  to  the  suppression  of  the  boundary 
layer,  viscosity  of  nanofluids  as  well  as  dispersion  of  the  nanoparti¬ 
cles.  Fig.  13  also  demonstrates  particle  volume  fraction  dependence 
of  Nusselt  number.  The  Nusselt  number  of  this  nanofluid  is  found 
to  increase  almost  linearly  with  the  particle  volume  fraction.  The 


nanofluid  behaves  more  like  a  fluid  than  a  conventional  solid 
(micrometer  or  millimeter)-fluid  mixture.  The  effects  of  several 
factors  such  as  gravity,  Brownian  force,  and  friction  force  between 
the  fluid  and  the  ultra-fine  particles  may  coexist  in  the  main  flow 
of  nanofluids. 


4.  Concluding  remarks 

A  comprehensive  review  on  three  major  cooling 
features— boiling,  droplet  spreading  and  convective  heat  transfer 
of  nanofluids  together  with  representative  results  from  own 
experimental  investigations  on  these  areas  are  presented  and 
analyzed  in  this  paper. 

From  the  review  of  available  experimental  results  for  boiling 
heat  transfer,  it  can  be  conferred  that  despite  of  contradictory 
and  inconsistent  data  on  boiling  heat  transfer,  there  is  undisputed 
substantial  increase  in  the  boiling  critical  heat  flux  of  nanofluids. 
However,  reported  data  are  still  limited  and  scattered  to  clearly 
understand  the  underlying  mechanisms  and  trend  of  boiling  heat 
transfer  performance  of  nanofluids.  In  addition,  only  a  couple  of 
efforts  are  made  on  flow  boiling  of  nanofluids.  Thus,  it  is  impera¬ 
tive  to  conduct  more  research  on  flow  boiling  of  nanofluids  under 
the  influence  of  various  factors  such  as  pressure,  mass  flux  and 
subcooling  besides  performing  more  systematic  experimental  and 
theoretical  investigations  on  pool  boiling  features  of  nanofluids. 
Without  providing  physical-chemistry  based  details  explanation, 
there  are  common  presumptions  such  as  deposition  of  nano¬ 
particles  or  tubes  on  heat  transfer  surface  and  surface  wettability 
used  as  reasons  for  the  observed  results  on  boiling  heat  transfer  of 
nanofluids. 

This  paper  also  presents  investigations  of  the  effects  of  surfac¬ 
tant  concentration  and  surface  tension  on  the  pool  boiling  heat 
transfer  of  carbon  nanotubes-nanofluid  and  results  showed  that 
large  enhancement  of  boiling  heat  flux  is  possible  and  would 
depend  on  the  concentration  of  the  surfactants.  The  effect  of  relax¬ 
ation  of  surface  tension  of  nanofluid  on  its  burnout  heat  flux  is  found 
to  be  significant. 

Although  spreading  of  liquid  droplet  plays  a  key  role  in  many 
industrial  processes  like  spray  cooling,  coating  and  ink-jet  printing, 
only  a  couple  of  research  efforts  have  been  made  to  study  droplet 
impingement  characteristics  of  nanofluids  on  solid  surfaces  under 
various  conditions.  Thus  more  studies  are  needed  on  dynamics  of 
non-boiling  and  boiling  nanofluids  droplets  impinging  on  solid  or 
liquid  surfaces. 

The  representative  results  and  the  review  of  the  findings 
from  the  literature  on  convective  heat  transfer  of  nanofluids 
demonstrated  that  nanofluids  exhibit  an  enhanced  heat  transfer 
coefficient  compared  to  its  base  fluid  and  it  increases  significantly 
with  increasing  concentration  of  nanoparticles  as  well  as  Reynolds 
number.  The  review  also  shows  a  considerable  chaos  and  random¬ 
ness  in  the  reported  data  on  convective  heat  transfer  coefficient 
from  various  research  groups.  Therefore,  more  careful  and  system¬ 
atic  investigations  on  nanofluids  preparation  and  the  convective 
heat  transfer  measurements  are  needed.  A  clear  understanding  of 
the  convective  heat  transfer  mechanisms  in  nanofluids  is  also  not 
yet  established. 

Researchers  working  on  nanofluids  have  shown  tremendous 
attention  to  the  static  thermal  properties,  particularly  thermal 
conductivity  of  nanofluids  rather  than  the  boiling  and  convective 
heat  transfer  features.  Although  the  applications  of  nanofluids  as 
advanced  cooling  medium  appears  to  be  promising,  the  advance¬ 
ment  toward  concrete  understanding  on  observed  properties  and 
features  of  nanofluids  as  well  as  their  development  for  commercial 
applications  remain  challenging  mainly  due  to  lack  of  agreement 
in  the  data  from  different  research  groups,  lack  of  understanding 
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of  the  mechanisms,  and  unsystematic  measurements  and  sample 
preparation.  Hence,  proper  sample  preparation  and  repeatable  and 
more  systematic  experimental  studies  on  measuring  any  properties 
of  nanofluids  are  worthwhile. 

In  renewable  energy  sector,  nanofluids  could  potentially  be  used 
to  enhance  heat  energy  storage  from  solar  collectors  and  to  increase 
the  energy  density.  Thus,  future  research  should  also  extend  to  such 
renewable  energy-based  applications  of  nanofluids. 

References 

[1]  Lide  DR.  CRC  handbook  of  chemistry  and  physics.  Taylor  and  Francis;  2007. 

[2]  Kaviany  M.  Principles  of  heat  transfer.  John  Wiley  &  Sons,  Inc.;  2002. 

[3  ]  Masuda  H,  Ebata  A,  Teramae  K,  Hishinuma  N.  Alteration  of  thermal  conductivity 
and  viscosity  of  liquid  by  dispersing  ultra-fine  particles  (dispersion  of  7-AI2O3, 
Si02,  and  Ti02  ultra-fine  particles).  Netsu  Bussei  1993;4:227-33. 

[4]  Slack  GA.  Thermal  conductivity  of  MgO,  AI2O3,  MgAl204,  and  Fe304  Crystals 
from  3  to  300  K.  Physical  Review  1962;126:427-41. 

[5]  Slack  GA.  Thermal  conductivity  of  pure  and  impure  silicon,  silicon  carbide  and 
diamond.  Journal  of  Applied  Physics  1964;35:3460-6. 

[6]  Kim  P,  Shi  L,  Majumdar  A,  McEuen  PL.  Thermal  transport  measurements  of 
individual  multiwalled  nanotubes.  Physical  Review  Letters  2001  ;87:215502- 
1-4. 

[7]  Yang  YM,  MaaJR.  Boiling  of  suspension  of  solid  particles  in  water.  International 
Journal  of  Heat  and  Mass  Transfer  1984;27:145-7. 

[8]  Choi  SUS.  Enhancing  thermal  conductivity  of  fluids  with  nanoparticles.  ASME 
FED  1995;231:99-105. 

[9]  Gass  V,  Van  der  Schoot  BH,  de  Rooij  NF.  Nanofluid  handling  by  micro-flow- 
sensor  based  on  drag  force  measurements.  In:  Proceedings  of  IEEE  conference 
on  micro  electro  mechanical  systems.  1993.  p.  167-72. 

[10]  Grimm  A.  Powdered  aluminum-containing  heat  transfer  fluids.  German  patent 
DE  4131516  Al;  1993. 

[11]  Lee  S,  Choi  SUS,  Li  S,  Eastman  JA.  Measuring  thermal  conductivity  of  fluids 
containing  oxide  nanoparticles.  Journal  of  Heat  Transfer  1999;121:280-9. 

[12]  Murshed  SMS,  Leong  KC,  Yang  C.  Enhanced  thermal  conductivity  of  Ti02  -water 
based  nanofluids.  International  Journal  of  Thermal  Sciences  2005 ;44: 
367-73. 

[13]  Murshed  SMS,  Leong  KC,  Yang  C.  Investigations  of  thermal  conductivity  and 
viscosity  of  nanofluids.  International  Journal  of  Thermal  Sciences  2008 ;47: 
560-8. 

[14]  Yu  W,  France  DM,  Routbort  JL,  Choi  SUS.  Review  and  comparison  of  nanofluid 
thermal  conductivity  and  heat  transfer  enhancements.  Heat  Transfer  Engineer¬ 
ing  2008;29:432-60. 

[15]  Murshed  SMS,  Leong  KC,  Yang  C.  Thermophysical  and  electrokinetic 
properties  of  nanofluids— a  critical  review.  Applied  Thermal  Engineering 
2008;28:2109-25. 

[16]  Murshed  SMS,  Leong  KC,  Yang  C.  Determination  of  the  effective  thermal  dif- 
fusivity  of  nanofluids  by  the  double  hot-wire  technique.  Journal  of  Physics  D: 
Applied  Physics  2006;39:5316-22. 

[17]  Buongiorno,  et  al.  A  benchmark  study  on  the  thermal  conductivity  of  nanoflu¬ 
ids.  Journal  of  Applied  Physics  2009;106:094312. 

[18]  Tong  LS,  Tang  YS.  Boiling  heat  transfer  and  two-phase  flow.  Taylor  &  Francis; 
1997. 

[19]  Pioro  IL,  Rohsenow  W,  Doerffer  SS.  Nucleate  pool-boiling  heat  transfer.  I. 
Review  of  parametric  effects  of  boiling  surface.  International  Journal  of  Heat 
and  Mass  Transfer  2004;47:5033-44. 

[20]  Pioro  IL,  Rohsenow  W,  Doerffer  SS.  Nucleate  pool-boiling  heat  transfer.  II. 
Assessment  of  prediction  methods.  International  Journal  of  Heat  and  Mass 
Transfer  2004;47:5045-57. 

[21]  Kang  MG.  Experimental  investigation  of  tube  length  effect  on  nucleate  pool 
boiling  heat  transfer.  Annals  of  Nuclear  Energy  1998;25:295-304. 

[22  ]  Kang  MG.  Effect  of  surface  roughness  on  pool  boiling  heat  transfer.  International 
Journal  of  Heat  and  Mass  Transfer  2000;43:4073-85. 

[23]  Kang  MG.  Effects  of  tube  inclination  on  pool  boiling  heat  transfer.  Nuclear 
Engineering  and  Design  2003;220:67-81. 

[24]  You  SM,  Kim  JH,  Kim  KM.  Effect  of  nanoparticles  on  critical  heat  flux  of  water 
in  pool  boiling  of  heat  transfer.  Applied  Physics  Letters  2003;83:3374-6. 

[25]  Milanova  D,  Kumar  R.  Heat  transfer  behavior  of  silica  nanoparticles  in  pool 
boiling  experiment.  Journal  of  Heat  Transfer  2008;130:042401-l-6. 

[26]  Das  SK,  Putra  N,  Roetzel  W.  Pool  boiling  characterization  of  nano-fluids.  Inter¬ 
national  Journal  of  Heat  and  Mass  Transfer  2003;46:851-62. 

[27]  Wen  D,  Ding  Y.  Experimental  investigation  into  the  pool  boiling  heat  trans¬ 
fer  of  aqueous  based  alumina  nanofluids.  Journal  of  Nanoparticle  Research 
2005;7:265-74. 

[28]  Witharana  S.  Boiling  of  refrigerants  on  enhanced  surfaces  and  boiling  of 
nanofluids.  PhD  thesis.  Royal  Institute  of  Technology;  2003. 

[29]  Prakash  NG,  Anoop  KB,  Das  SK.  Mechanism  of  enhancement/deterioration  of 
boiling  heat  transfer  using  stable  nanoparticles  suspensions  over  vertical  tubes. 
Journal  of  Applied  Physics  2007;102:074317-l-7. 

[30]  Prakash  NG,  Anoop  KB,  Sateesh  G,  Das  SK.  Effect  of  surface  orientation  on 
pool  boiling  heat  transfer  of  nanoparticle  suspensions.  International  Journal 
of  Multiphase  Flow  2008;34:145-60. 


[31]  Soltani  S,  Etemad  SG,  Thibault  J.  Pool  boiling  heat  transfer  performance  of 
Newtonian  nanofluids.  Heat  and  Mass  Transfer  2009;45:1555-60. 

[32]  Soltani  S,  Etemad  SG,  Thibault  J.  Pool  boiling  heat  transfer  performance  of  New¬ 
tonian  nanofluids.  International  Communications  of  Heat  and  Mass  Transfer 
2010;37:29-33. 

[33]  Truong  B,  Hu  LW,  Buongiorno  J,  McKrell  T.  Modification  of  sandblasted  plate 
heaters  using  nanofluids  to  enhance  pool  boiling  critical  heat  flux.  International 
Journal  of  Heat  and  Mass  Transfer  2010;53:85-94. 

[34]  Kim  SJ,  McKrell  T,  Buongiorno  J,  Hu  LW.  Enhancement  of  flow  boiling  crit¬ 
ical  heat  flux  (CHF)  in  alumina/water  nanofluids.  Advanced  Science  Letters 
2009;2:100-2. 

[35]  Henderson  K,  Park  YG,  Liu  L,  Jacobi  AM.  Flow-boiling  heat  transfer  of  R-134a- 
based  nanofluids  in  a  horizontal  tube.  International  Journal  of  Heat  and  Mass 
Transfer  2010;53:944-51. 

[36]  Vassallo  P,  Kumar  R,  D’Amico  S.  Pool  boiling  heat  transfer  experiments  in 
silica-water  nano-fluids.  International  Journal  of  Heat  and  Mass  Transfer 
2004;47:407-11. 

[37]  Bang  IC,  Chang  SH.  Boiling  heat  transfer  performance  and  phenomena  of 
Al203-water  nano-fluids  from  a  plain  surface  in  a  pool.  International  Journal 
of  Heat  and  Mass  Transfer  2005;48:2407-19. 

[38]  Kim  H,  Kim  J,  Kim  M.  Experimental  study  on  CHF  characteristics  of  water-Ti02 
nano-fluids.  Nuclear  Engineering  Technology  2006;39:61-8. 

[39]  Kim  SL,  Bang  IC,  Buongiorno  J,  Hu  LW.  Surface  wettability  change  during  pool 
boiling  of  nanofluids  and  its  effect  on  critical  heat  flux.  International  Journal  of 
Heat  and  Mass  Transfer  2007;50:4105-16. 

[40]  Jackson  J.  Investigation  into  the  pool-boiling  characteristics  of  gold  nanofluids. 
MS  thesis.  University  of  Missouri-Columbia;  2007. 

[41  ]  Chopkar  M,  Das  AK,  Manna  I,  Das  PK.  Pool  boiling  heat  transfer  characteristics 
of  Zr02-water  nanofluids  from  a  flat  surface  in  a  pool.  Heat  and  Mass  Transfer 
2008;44:999-1004. 

[42]  Liu  ZH,  Liao  L.  Sorption  and  agglutination  phenomenon  of  nanofluids  on  a  plain 
heating  surface  during  pool  boiling.  International  Journal  of  Heat  and  Mass 
Transfer  2008;51 :2593-602. 

[43]  Lv  LC,  Liu  ZH.  Boiling  characteristics  in  small  vertical  tubes  with  closed  bot¬ 
tom  for  nanofluids  and  nanoparticles-suspensions.  Heat  and  Mass  Transfer 
2008;45:1-9. 

[44]  Kathiravan  R,  Kumar  R,  Gupta  A,  Chandra  R,  Jain  PK.  Pool  boiling  characteristics 
of  carbon  nanotube  based  nanofluids  over  a  horizontal  tube.  Journal  of  Thermal 
Science  and  Engineering  Applications  2009;  1 :022001 -1-7. 

[45]  Rohsenow  WM.  A  method  of  correlating  heat  transfer  data  for  surface  boiling 
of  liquids.  Transactions  of  the  ASME  1952;74:969-76. 

[46]  Zuber  N.  Hydrodynamic  aspects  of  boiling  heat  transfer.  Physics  and  Mathe¬ 
matics,  AEC  Report  No.  AECU-4439;  1959. 

[47]  Pioro  IL.  Experimental  evaluation  of  constants  for  the  Rohsenow  pool  boil¬ 
ing  correlation.  International  Journal  of  Heat  and  Mass  Transfer  1999;42: 
2003-13. 

[48]  Murshed  SMS,  Milanova  D,  Kumar  R.  An  experimental  study  of  surface  tension- 
dependent  pool  boiling  characteristics  of  carbon  nanotubes-nanofluids. 
In:  Proceedings  of  7th  international  ASME  conference  on  nanochannels, 
microchannels  and  minichannels.  2009. 

[49]  Milanova  D,  Kumar  R.  Functionalized  single  walled  and  double  walled  carbon 
nanotubes  for  thermal  enhancement.  In:  Proceedings  of  the  ASME  international 
mechanical  engineering  congress  and  exposition.  2007. 

[50]  Kim  H,  Kim  J,  Kim  MH.  Effect  of  nanoparticles  on  CHF  enhancement  in 
pool  boiling  of  nano-fluids.  International  Journal  of  Heat  and  Mass  Transfer 
2006;49:5070-4. 

[51]  Wasan  DT,  Nikolov  AD.  Spreading  of  nanofluids  on  solids.  Nature 
2003;423:156-9. 

[52]  Duursma  G,  Sefiane  K,  Kennedy  A.  Experimental  studies  of  nanofluid  droplets 
in  spray  cooling.  Heat  Transfer  Engineering  2009;30:1108-20. 

[53]  Shen  J,  Liburdy  JA,  Pence  DV,  Narayanan  V.  Droplet  impingement  dynamics: 
effect  of  surface  temperature  during  boiling  and  non-boiling  conditions.  Journal 
of  Physics:  Condensed  Matter  2009;21 :464133-1-14. 

[54]  Murshed  SMS,  Nieto  de  Castro  CA.  Spreading  characteristics  of  nanofluid 
droplets  impacting  onto  a  solid  surface.  Journal  of  Nanoscience  and  Nanotech¬ 
nology  2011;11:1-7. 

[55]  Pak  BC,  Cho  YI.  Hydrodynamic  and  heat  transfer  study  of  dispersed  flu¬ 
ids  with  submicron  metallic  oxide  particles.  Experimental  Heat  Transfer 
1998;11:151-70. 

[56]  Eastman  JA,  Choi  SUS,  Li  S,  Soyez  G,  Thompson  LJ,  Dimelfi  RJ.  Novel 
thermal  properties  of  nanostructured  materials.  Materials  Science  Forum 
1999:312-314:629-34. 

[57]  Xuan  Y,  Li  Q,  Investigation  on  convective  heat  transfer  and  flow  features  of 
nanofluids.  Journal  of  Heat  Transfer  2003;125:151-5. 

[58]  Wen  D,  Ding  Y.  Experimental  investigation  into  convective  heat  transfer  of 
nanofluids  at  the  entrance  region  under  laminar  flow  conditions.  International 
Journal  of  Heat  and  Mass  Transfer  2004;47:5181-8. 

[59]  Heris  SZ,  Etemad  SG,  Esfahany  MS.  Experimental  investigation  of  oxide  nanoflu¬ 
ids  under  laminar  flow  convective  heat  transfer.  International  Communications 
of  Heat  and  Mass  Transfer  2006;33:529-35. 

[60]  Duangthongsuk  W,  Wongwises  S.  Heat  transfer  enhancement  and  pressure 
drop  characteristics  of Ti02  -water  nanofluid  in  a  double-tube  counter  flow  heat 
exchanger.  International  Journal  of  Heat  and  Mass  Transfer  2009;52:2059-67. 

[61  ]  Faulkner  D,  Rector  DR,  Davison  JJ,  Shekarriz  R.  Enhanced  heat  transfer  through 
the  use  of  nanofluids  in  forced  convection.  In:  Proceedings  of  ASME  heat  trans¬ 
fer  div.  2004.  p.  219-24. 


2354 


S.M.S.  Murshed  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  15  (201  \ )  2342-2354 


[62]  Jung  JY,  Oh  HS,  Kwak  HY.  Forced  convective  heat  transfer  of  nanofluids  in 
microchannels.  In:  Proceedings  of  the  ASME  international  mechanical  engi¬ 
neering  congress  and  exposition.  2006. 

[63]  Ding  Y,  Alias  H,  Wen  D,  Williams  AR.  Heat  transfer  of  aqueous  suspen¬ 
sions  of  carbon  nanotubes.  International  Journal  of  Heat  and  Mass  Transfer 
2006;49:240-50. 

[64]  Yang  Y,  Zhang  ZG,  Grulke  EA,  Anderson  WB,  Wu  G.  Heat  transfer  properties 
of  nanoparticle-in-fluid  dispersions  (nanofluids)  in  laminar  flow.  International 
Journal  of  Heat  and  Mass  Transfer  2005;48:1107-16. 

[65]  Lai  WY,  Duculescu  B,  Phelan  PE,  Prasher  RS.  Convective  heat  transfer  with 
nanofluids  in  a  single  1 .02-mm  tube.  In:  Proceedings  of  the  ASME  international 
mechanical  engineering  congress  and  exposition.  2006. 

[66]  Williams  W,  Buongiorno  J,  Hu  LW.  Experimental  investigation  of  turbulent  con¬ 
vective  heat  transfer  and  pressure  loss  of  alumina/water  and  zirconia/water 
nanoparticle  colloids  (nanofluids)  in  horizontal  tubes.  Journal  of  Heat  Transfer 
2008;130:1-7. 

[67]  Hwang  KS,  Jang  SP,  Choi  SUS.  Flow  and  convective  heat  transfer  characteris¬ 
tics  of  water-based  AI2O3  nanofluids  in  fully  developed  laminar  flow  regime. 
International  Journal  of  Heat  and  Mass  Transfer  2009;52:193-9. 

[68]  Xie  H,  Li  Y,  Yu  W.  Intriguingly  high  convective  heat  transfer  enhancement  of 
nanofluid  coolants  in  laminar  flows.  Physics  Letters  A  2010;374:2566-8. 

[69]  Amrollahi  A,  Rashidi  AM,  Loth  R,  Meibodi  ME,  Kashefi  K.  Convection  heat  trans¬ 
fer  of  functionalized  MWNT  in  aqueous  fluids  in  laminar  and  turbulent  flow  at 
the  entrance  region.  International  Communications  in  Heat  and  Mass  Transfer 
2010;37:717-23. 

[70]  Kulkarni  DP,  Namburu  PK,  Bargar  HE,  Das  DK.  Convective  heat  transfer  and  fluid 
dynamic  characteristics  of  Si02 -ethylene  glycol/water  nanofluid.  Heat  Transfer 
Engineering  2008;29:1027-35. 

[71  ]  Daungthongsuk  W,  Wongwises  S.  A  critical  review  of  convective  heat  transfer 
of  nanofluids.  Renewable  and  Sustainable  Energy  Reviews  2007;11:797-817. 

[72]  Kaka£  S,  Pramuanjaroenkij  A.  Review  of  convective  heat  transfer  enhance¬ 
ment  with  nanofluids.  International  Journal  of  Heat  and  Mass  Transfer 
2009;52:3187-96. 

[73]  Putra  N,  Roetzel  W,  Das  SK.  Natural  convection  of  nanofluids.  Heat  and  Mass 
Transfer  2003;39:775-84. 

[74]  Khanafer  K,  Vafai  K,  Lightstone  M.  Buoyancy-driven  heat  transfer  enhancement 
in  a  two-dimensional  enclosure  utilizing  nanofluids.  International  Journal  of 
Heat  and  Mass  Transfer  2003;46:3639-53. 

[75]  Wen  D,  Ding  Y.  Formulation  of  nanofluids  for  natural  convective  heat  transfer 
applications.  International  Journal  of  Heat  and  Fluid  Flow  2005;26:855-64. 

[76]  Buongiorno  J.  Convective  transport  in  nanofluids.  Journal  of  Heat  Transfer 
2006;128:240-50. 

[77 ]  Bejan  A.  Convection  heat  transfer.  John  Wily  &  Sons,  Inc.;  2004. 


[78]  Dittus  FW,  Boelter  LMK.  Heat  transfer  in  automobile  radiators  of  the  tubular 
type.  University  of  California  Publications  on  Engineering  1930;2:443-61. 

[79]  Mai'ga  EEB,  Palm  SJ,  Nguyen  CT,  Roy  G,  Galanis  N.  Heat  transfer  enhancement 
by  using  nanofluids  in  forced  convection  flows.  International  Journal  of  Heat 
and  Fluid  Flow  2005;26:530-46. 

[80]  Roy  G,  Nguyen  CT,  Lajoie  PR.  Numerical  investigation  of  laminar  flow  and  heat 
transfer  in  a  radial  flow  cooling  system  with  the  use  of  nanofluids.  Superlattices 
and  Microstructures  2004;35:497-511. 

[81  ]  Xuan  Y,  Yu  K,  Li  Q,  Investigation  on  flow  and  heat  transfer  of  nanofluids  by  the 
thermal  Lattice  Boltzmann  model.  Progress  in  Computational  Fluid  Dynamics 
2005;5:13-9. 

[82]  Palm  SJ,  Roy  G,  Nguyen  CT.  Heat  transfer  enhancement  with  the  use  of 
nanofluids  in  radial  flow  cooling  systems  considering  temperature-dependent 
properties.  Applied  Thermal  Engineering  2006;26:2209-18. 

[83]  Jou  RY,  Tzeng  SC.  Numerical  research  of  nature  convective  heat  transfer 
enhancement  filled  with  nanofluids  in  rectangular  enclosures.  International 
Communications  in  Heat  and  Mass  Transfer  2006;33:727-36. 

[84]  Jang  SP,  Choi  SUS.  Cooling  performance  of  a  microchannel  heat  sink  with 
nanofluids.  Applied  Thermal  Engineering  2006;26:2457-63. 

[85]  Heris  SZ,  Esfahany  MN,  Etemad  G.  Numerical  investigation  of  nanofluid  laminar 
convective  heat  transfer  through  a  circular  tube.  Numerical  Heat  Transfer  Part 
A:  Applications  2007;52:1043-58. 

[86]  Mirmasoumi  S,  Behzadmehr  A.  Effect  of  nanoparticles  mean  diameter  on  mixed 
convection  heat  transfer  of  a  nanofluid  in  a  horizontal  tube.  International  Jour¬ 
nal  of  Heat  and  Fluid  Flow  2008;29:557-66. 

[87]  Ho  CJ,  Chen  MW,  Li  ZW.  Numerical  simulation  of  natural  convection  of 
nanofluid  in  a  square  enclosure:  effects  due  to  uncertainties  of  viscosity 
and  thermal  conductivity.  International  Journal  of  Heat  and  Mass  Transfer 
2008;51:4506-16. 

[88]  Oztop  HF,  Abu-Nada  E.  Numerical  study  of  natural  convection  in  partially 
heated  rectangular  enclosures  filled  with  nanofluids.  International  Journal  of 
Heat  and  Fluid  Flow  2008;29:1326-36. 

[89]  Vajjha  RS,  Das  DK,  Namburu  PK.  Numerical  study  of  fluid  dynamic  and  heat 
transfer  performance  of  Al2  O3  and  CuO  nanofluids  in  the  flat  tubes  of  a  radiator. 
International  Journal  of  Heat  and  Fluid  Flow  2010;31:613-21. 

[90]  Murshed  SMS,  Leong  KC,  Yang  C,  Nguyen  NT.  Convective  heat  transfer 
characteristics  of  aqueous  Ti02  nanofluids  under  laminar  flow  conditions. 
International  Journal  of  Nanoscience  2008;7:325-31. 

[91]  Hamilton  RL,  Crosser  OK.  Thermal  conductivity  of  heterogeneous  two 
component  systems.  Industrial  and  Engineering  Chemistry  Fundamentals 
1962;1:187-91. 

[92]  Batchelor  GK.  The  effect  of  Brownian  motion  on  the  bulk  stress  in  a  suspension 
of  spherical  particles.  Journal  of  Fluid  Mechanics  1977;83:97-117. 


